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Abstract. In this note, we are concerned with the asymptotic approximation of a class 
of double integrals which can be represented as an angular spectrum superposition. These 
double integrals typically appear in electromagnetic scattering problems. Based on the 
synthetic manipulation of the method of steepest descent path, approximate expressions 
of the double integrals are derived in terms of the leading term of the contribution to the 
asymptotic expansions. 

1 Introduction 

The methods of steepest descent path (SDP) and stationary phase are commonly useful 
for the asymptotic evaluation of the complicated integrals in electromagnetics particularly 
in the far zoneLang-LM. For the one-dimensional case, it seems no restrictions need 
to be enforced on the integrands to apply these classic asymptotic methods, although 
sometimes the construction of the SDP is difficult as also is the establishment of the 
uniform asymptotic expansions in the neighborhoods of the singularities of the integrands 
close to the saddle points Wong, Van. But for the asymptotic approximation of double 
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integrals 

liX) = JJ (pix,y)exp[i\h{x,y)]dxdy , X^oo, (1) 

D 

the stationary pliasc method is gcneraUy used provided tliat the integrands liave certain 
restricted forms, i.e., the phase functions h{x, y) are real-valued BW, Wong, BMM, Jonel. 
When the complex-valued phase functions are concerned, asymptotic expansions of the 
double integrals are still possible to be obtained if the integrands are smooth sufficiently 
and vanish on the boundary of the integration domain D. In this case, the dominant 
contribution to the asymptotic expansions may come from the critical points of the first 
kind. However, if the phase function is complex-valued in general but not differentiable 
everywhere in the domain, some additional work may be necessary to account for the con- 
tribution of the boundary stationary points to the asymptotic expansions of the double 
integrals Wongl-Jone2. These boundary stationary points are certainly located in the 
vicinity of regions wherein the phase function is not differentiable, i.e., branch points and 
singularities. Therefore a curve of stationary points may be formed on the boundaries of 
the sub-domains Dj C D which are disjoined with each other by the "deleted" neighbor- 
hoods of nondifferentiation. In addition, because the phase functions now are complex- 
valued, approximation of the Laplace-type integrals may need to be manipulated on some 
sub-domains to accomplish the classic asymptotic evaluation of these double integrals 
Wong,^Jone2 

In this note we present a simpler way to establish the asymptotic approximation of a 
class of double integrals which are derived from an angular spectrum representation of 

the field-related entities. After departing from the SDP method operated in a synthetic 

^The boundary stationary points are typically the critical points of the second kind. If there is a curve 
of stationary points formed on the boundary of the domain D, we know that the leading term of the 
contribution to the asymptotic expansion of the double integral in (1) is 0(A^^) as A ^ oo. Similarly, 
for the asymptotic expansion of the Laplace type integral J(A) = // j/) exp [— A /i(a;, j/)] da;dj/, if the 

D 

stationary points are the critical points of the first kind, then the leading term of /(A) is also 0(A~^) as 
A — > +00. For more details, see Wong,^Jone2. 
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fashion, the approximation of these double integrals is achieved by using the contour 
deformation to confine the integration on a small domain defined by two local SDPs. The 
saddle points along each SDP are employed to construct the truncated Taylor expansion 
of the complex-valued phase function which is differentiable on the mapped domain. The 
double integral is then asymptotically evaluated by the lowest order in the far-zone limit. 
A condition for the asymptotic approximation to be valid is also briefly commented. 



2 Theory 

To begin with let us suppose the vector r = (x, y, z) deflned on the upper half space 
= {z > 0} and a double integral G(r) generally having the form 

G{r) = J J f(kx, ky) exp [i{kxX + ky-y + k^z)] dk^ dky , (2) 

where kz = k^ — k"^ — k"^ is defined on the top Riemann surface such that lm.{kz) > for 
k'^+ky G M. Generally speaking, G(r) represents a class of functions that have the physical 
importance in reducing a field-related entity into an angular spectrum superposition. It 
is reasonably assumed that |G(r)| < oo is well defined by the integrand belonging to 
L^(M?). In most cases, we may also assume that f{kx, ky) in a21 is independent of r and 
has no singularities on the complex domain f2 = S)x]R'^U]R'^x2), where 2) is a complex 
neighborhood of the real axis, and M"^ = C \ M. 

According to Pubini's theorem, G{r) in a21 could be rewritten as 

/oo /•oo 
e-xp{ikxx) dkx f{kx,ky)exp[i{kyy + kzz)]dky. (3) 
-oo J —oo 



-oo 
oo 



By denoting g{kx) = J f{kx, ky) exp [i{kyy + kzz)] dky, we know that g{kx) G L''^(M) al- 

— oo 

most everywhere. Indeed, g{kx) may have isolated singularities which depend on f{kx, ky). 
However, because f{kx, ky) is assumed to be nonsingular in fi, it is clear that the isolated 
singularities oi g{kx)i if exist, are only real— valued. Conventionally g{kx 

) can be expanded 



in the form g{kx) = {kx — Sn)^" ^ in the vicinity of the singular point Sn G M. Since 

n=0 

g{kx) G i^^(lK), it seems that those singularities {s„} are not likely to be poles, so that 
jn > 0. On the other hand, we can reform g{kx) as the multiplication of two parts 

9{kx) = f{kx) expiikpxPx) , (4) 

where px = Vr"^ — x^, r = |r| and kpx = ^/k^ — k^ with lm.{kpx) > for real kx- By 
substituting a23 into a22 and noting that }(kx) has no poles, it is always feasible to reduce 
G(r) in a22 into a contour integral along the SDP of the phase function. Therefore, we 
obtain synthetically the first asymptotic approximation 



/oo 
f{kx) exp [i{kxx + kpxPx)] ^kx 
-oo 

= j Rkx) exp [ikorcosioLx - il^x)] 



SDPi 

~ j Rkx) exp [ikorcos{ax - ipx)] dkx , (5) 

as kor — oo. In a24, kx = kocosax, ipx = cos~^ (^) and SDPi : (—00,00) 1— > C is deter- 
mined by the parametrization of ax (or kx) through the equation Re [cos(Q;a; — -^a;) — 1] = 
0. SBF]"*^ denotes the local path of the SDPi that passes through the saddle point 

kxs = ko COsV'a;. 

On account of a22-a24, by using Fubini's theorem again we obtain 

/oo p 
exp{ikyy)dky / f{kx,ky) e-xp[i{kxX + k^z)] dkx 
-00 J 

sDP'r 

/oo ~ 
f{ky) exp [i{kyy + kpyPy)] dky , (6) 
-00 

where the function f{ky) is defined by the equation 

J f{kx, ky) exp [i{kxX + k^z)] dkx = f{ky) ey:p{ikpyPy) . (7) 
SDP'r 



Here py = \lr'^ — and kpy = y^o ~ "^ith \m.{kpy) > for real ky. Actually the 
definition of f{ky) in a26 is very similar to that of f{kx) in a23. Following the same 
approach as illustrated in a22-a24, we can obtain the second asymptotic approximation 
of G{r) as follows: 

G(r) J f {ky) exp [ikorcos{ay — ipy)] dky 
SDP2 

f{ky) exp [ikorcos{ay - ipy)] dky . (8) 
SDP2°" 

In a27, ky = kocosay, i/jy = cos~^ (^) and SDP2 : (—00,00) 1— >■ C is determined by the 
parametrization of ay (or ky) through the equation Re [cos(a2^ — ijjy) — 1] =0. SDPg"^ 
denotes the local path of the SDP2 that passes through the saddle point kyg = ko cosipy. 

It seems that G(r) can be asymptotically approximated by either a24 or a27. Clearly, 
the asymptotic expression of G(r) is strongly dependent on the properties of f{kx) and 
f{ky) in the fact that the leading contribution to G(r) would be significantly influenced 
or even dominated by the singularities of f{kx) (or f{ky)) that are close to the saddle 
point kxs (or kyg) Wong, 3Vong2. On the other hand, the synthetic equations a23 and 
a26 imply that both f{kx) and f{ky) have been modulated by the larger parameter kor. 
Therefore, the asymptotic expansion of G{r) cannot be obtained in closed-form, unless 
the asymptotic expansion of f{kx) (or f{ky) ) is achieved beforehand. 

However, the asymptotic approximation of G(r) can still be accomplished by restoring 
G(r) with a double integral representation. In fact, due to a25-a27, the approximation of 
G(r) is now reduced to a double integral 

G{y) Ri j j f{kx, ky) exp [i{kxX + kyy + k^z)] dkx dky , (9) 

D' 

where the domain C Z) D' := SDP[°'' x SDP^"^ is located in a neighborhood of the critical 
point {kxs,kys). The local pathes of the SDPs can be represented parametrically by the 
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real-valued ^ and r] as 

kxiO = kxs + - i)^ , (10) 

ky{ii) = kys + kzs{l - i)r] , (11) 

where kzs = ^Jk^ — k%g — k^g Fels,enl. Therefore D' can be viewed a continuous map from 
the local domain ^ x r/ C by means of a29 and a210. 

Let U{(,r]) = {kor)~^{kxX + kyy + kzz) — 1 be defined on D'; then G(r) in a28 can be 
written in the form 

G(r) « exp{ikor) J J f{kx,ky)exp[ikQrU{(,,r])] d^dr/. 

D' 

~ fikxs,kys)exp{ikor) J J exp[ikQrU{$,,ri)] d^d??, (12) 

by utilizing the approximation f{kx, ky) f{kxs, kyg) on D' . Indeed, to be more accurate, 
we need to express f{kx,ky) in term of the Taylor expansion by assuming f{kx,ky) is 
analytic in a neighborhood of the critical point. There will be no additional difficulty 
in working out the asymptotic approximation of G{r) in a211 when f{kx,ky) can be 



represented by its Taylor expansion. Furthermore, if U{^, rj) oi kz = y kQ — k'^ — ky is 
differentiable on the domain D', then we can approximate U((,,ri) by its truncated Taylor 
expansion on the critical point {^,r]) = (0,0). In fact, kz = -\/fef is holomorphic on the 
complex domain C D D2 := {k^} continuously mapped from D' as shown in Figure 1, 
simply because D2 is located on a Reimann surface that neither crosses the branch cut nor 
encloses the branch point (0,0). Therefore U{^^vi) is analytic on D' in accordance with 
the chain rule of differentiation. Then ?7(^, 77) is approximated by the Taylor expansion 
at (0, 0) truncated to the second order as: 

U{i,r^)^UT{i,il) = i{ae + bif-2cir^), i^,v) & D' , (13) 

by recalling ^ = ^ = at (0, 0). Here 

a = l-4, b=l-^, c=^, (14) 
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are real-valued constants. 

After substituting a212 into a211, we then derive the asymptotic approximation ( kor 

oo) 

G(r) f{kxs,kys)exp{ikor) JJ exp [-kor{a('^ + brj"^ - 2c^ri)] d^dr/ 

D' 

oo oo 

^ f{kxs,kys)exp{ikor) J J exp [-kor{a$,^ + bri^ - 2c^ri)] d^dri , (15) 

— oo — oo 

which is the leading term of the factor {kor)~^ of the contribution to the asymptotic 
expansion of G(r). In fact, the integral on the right side of a214 can be easily calculated 
by eliminating the corss-product term ^t] through a linear transformation. 

It is worthy to note that the approximation of U{^,r)) by its truncated Taylor ex- 
pansion a212 can be reasonably achieved only under the condition that ^ = ^ > is 
not too small. In fact the Taylor expansion of U{S,, rj) at (0,0) takes the form ?7(^, rj) = 

oo 

E ui3^l^\='^ [l + Pi3{Q\ where = and /? = {Pi, (32) is a multiindex such 

\0\=2 

that \P\ = Pi + 02- If the coefficient up is kept as a factor of 0(1), then the function 
PffiO ~ 0(|^|/6')l^l-2 could be achieved at some P for any \P\ > 2. After considering 
the approximation of a214, ~ 0{l/kor) G D' is at least guaranteed. Therefore, 
9 > O j is needed for an appropriate approximation of U{^,r]) by a212 on the 

domain D' . In other words, the asymptotic approximation of G(r) by a214 becomes 
reasonable only when 9 > 9q is necessarily satisfied, where = ■\J j^- 

3 Conclusion 

In this note, we are concerned with the asymptotic approximation of a class of double 
integrals which can be represented as an angular spectrum superposition. These double 
integrals are typically useful for the calculation of field-related problems in physics and 
electromagnetics. A simpler way is found to make the asymptotic approximation feasible 



and reasonable by applying the steepest descent path method on the double integral in 
a synthetic fashion. After reduction to a local integration on a small domain in which 
the complex-valued phase function is differentiable, the asymptotic approximation of the 
double integral is finally derived in terms of the leading term of the contribution to the 
asymptotic expansion. The validity of the asymptotic approximation is only under the 
condition that - > 9o is satisfied, where = \/ tt- > 0. 
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Figure Captions 



Figure 1. Schematic of the continuous map k'^ : D'{^,r)) D2{kl). 
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Figure 1: Schematic of the continuous map kj, : D'{^,ri) ^ D2{kj.). 
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